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Abstract 

Exact static non-degenerate stripe and checkerboard ground states are obtained in a two di- 
mensional generalized periodic Anderson model, for a broad concentration range below quarter 
filling. The random droplet states also present in the degenerate ground state, are eliminated by 
extending the Hamiltonian with terms of different physical origin such as dimerization, periodic 
charge displacements, density waves or distortion lines. 
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I. INTRODUCTION 



Quasi-one-dimensional intrinsic inhomogeneities called stripes 1 represent of the most chal- 
lenging problems in understanding self-organized structures ever since the discovery of high 
T c cuprates^. Such textures have been observed in other than superconducting cuprates, 
most notably in manganites 3 , nickelates 4 and rare-earth compounds 5 . Most recently it has 
been shown that all cuprates exhibit intrinsic inhomogeneities in some form or another: 
checkerboard structures have been identified in lightly hole-doped copper oxides 6 , while the 
electron-doped cuprates show droplet formations^- These textures are so predominant in 
high T c s, that it has been even suggested*^! as the best candidate for the electron pairing 
mechanism in these materials. 

After the existence of stripes have been established in a broad spectrum of materials, it be- 
comes clear that they emerge due to a microscopic incompatibility between two phases ^ 6 i 9 i 10 . 
However, this interplay not always leads to stripe formations, but could also lead to droplets 
(or blobs) 7,8 as well. Macroscopic inhomogeneity cannot occur as long as long-range forces 
are present. But when such forces are absent, depending on doping, local inhomogeneities 
appear in form of stripes or other form of clustering, e.g., droplets&ii. 

Even thought stripes have been known to exist much before high T c has been discovered^, 
their theoretical understanding lacks rigorous description. Hereafter we want to fill this gap. 
Hence we present an exact solution of a generalized periodic Anderson model (PAM) with 
ground states exhibiting intrinsic inhomogeneities of the stripe, checkerboard and droplets 
type. Using a two band model as a starting point of our analyses it even renders the obtained 
results to be more generally applicable, since real materials are mostly of multi-band type. 
These models are usually addressed by projecting the multi-band structure into a few-band 
picture^ 3 -, which we stop for mathematical convenience, at a two band level. However, our 
study is not a simple two-band model but rather an extended PAM which contains the 
added feature of strong correlation effects originating from the on-site Coulomb repulsion 
U present in the correlated (/) band, an arbitrary Ua > extra Hubbard interaction in the 
free (d) band, leaves our results unchanged. 

Since stripes and checkerboards are observed in a broad spectrum of materials, we are 
primarily focusing on ground states which exhibit these inhomogeneities, being less interested 
in the properties of the homogeneous phases in which they appear. 
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Our exact results can be summarized as follows. Below quarter filling, two stripe ground 
states emerge. One (/) with insulating and paramagnetic stripes, while the second one 
(II) with itinerant and ferromagnetic stripe lines. This second solution allows checkerboard 
structures as well. In both cases, i.e., I and //, the inter-stripe line regions contain empty 
sites, hence are insulating. The obtained ground states are in general degenerate. The 
degeneracy is provided by a random blob structure corresponding to the same energy, the 
blobs (random shape clusters) possessing the same properties as the stripe lines in cases 
I and 77. The degeneracy of the ground state can however be lifted, the resulting non- 
degenerate ground state remaining of pure stripe or checkerboard character. The lifting 
factor we find may have different physical origin as distortion lines, dimerization or periodic 
charge displacement (density waves). The obtained stripe formation processes are generic 
and are less sensitive to the properties of the homogenous phases present at quarter filling. 
We further note that marginal to the stripe problem, but interesting for the PAM itself, we 
were able to prove rigorously that the studied 2D PAM is ferromagnetic at quarter filling in 
a restricted region of its parameter space. Similar result has been recently reported for 3D 
PAM as weltii. 

The remaining part of the paper is structured as follows. Sec. II presents the exact trans- 
formation of the Hamiltonian in a positive semidefinite form, Sec. III. describes the obtained 
ground state solutions, Sec. IV. presents a discussion of the obtained results, and Sec.V. 
concluding the paper, closes the presentation. 



II. THE TRANSFORMATION OF THE HAMILTONIAN 

Our starting Hamiltonian is thus, H = H + U, written for a free (d) and a correlated 
(/) band, where we allow for hopping in both bands: 



£ t P Aji+r,a + V r (dtJ i+ria + 



P=dJ 



flJ i+r ,„) + H.c] + V Q (d\Ji,a + H.C.) + E f h{ a } . (1) 



To keep the description as general as possible, the above Hamiltonian is defined on a 2D 
Bravais lattice with unit cell X and primitive vectors (xi, X2). In Eq. V r , Vq, and 

Ef, represent the hopping amplitudes for p = d, f electrons, the non-local and on-site 
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hybridization, and the local / electron energy, respectively. The coordinate r, with possible 
values xi,X2,X2 ± xi, is allowed to extend to all sites of I. Within the unit cell X\ defined 
at i, the lattice sites rj. = i + r a/3 , with r a p = axi + /3x 2 , a, /3 = 0, 1, can be numbered by 
n(a,P) = 1 + a + 3/3 — 2a/3 without reference to Xi, see Fig. [T]a). The correlated / band 
has in addition U = U Ya ^{t^{j.> acting on it with < U < oo on-site Hubbard repulsion. 

It is known that exact solutions exist mostly in ID systems and in higher dimensions 
is almost impossible to find rigorous results. Thus, to find the exact ground state of H 
in 2D we use a method which is based on transforming H into a representation which is 
positive semidefinite. Quantum mechanics tells us that the minimum possible eigenvalue 
of a positive semidefinite operator, e.g., O, is zero. Hence the ground state, e.g., \^f g ) of 
O can be constructed from 0\^ g ) = 0. The fact that ground states containing intrinsic 
inhomogeneities can be obtained in this manner for a ,,non-integrable"— model as 2D-PAM 
does not come as a surprise, since the applied procedure works even in disordered case^. We 
further note that the method is described in details in Ref. Q and has been previously used 
to solve generalized PAM type models at 3/4 filling in 2D^£ and even in 3D X ^. In the present 
case the transformation is performed in 2D at a lower filling region, which is described in 
2D for the first time in this paper. 

The transformation proceeds in the following way: first we transform exactly the starting 
Hamiltonian in a positive semidefinite form. This is accomplished with the use of the 
operators 

A> = X! X! a n(a,l3),pPl+r a p,ai (2) 
a,0=O,lp=d,f 

which are linear combination of the original fermionic operators acting on the corners of an 
elementary plaquette , see Fig [I] a). It can be easily seen that P = J2i,a M. a ^-i,a contains 
exactly the same operators from Eq. Hence, properly choosing the coefficients a njP , 
n = n(a, 0), H from Eq. (JTJ can be written as H = P + E g , where E g is a constant. The 
proper mapping is 

Im 

£ a mj , p a* n . tp , = 7*f(t?, if, V r , E f ) , (3) 
i=i 

where for the r = anxi + 0:2X2 values allowed by (1), and p,p' = d, f, T£f is given by 

TPf = (1 - M[<W? + (1 - <W)K] + 6 v>0 [5 P;P ,(5 P)d K + 8 pJ {E f + K)) + (1 - 6 pj/ )V ), 
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nj = j5u,o + (j + \oci - a 2 \)5 Ut2 + (Sj )2 + (4|ai| + 2\a 2 \)5j ; i)8 u>1 , 
mj - nj = 26 u>2 + S u>1 [7 - 4 J + (6J - 10)|«i| + (67 - 8)|a 2 |]. 

21 M = 8 - 5u + u 2 , i/=(\a 1 \ + \a 2 \). (4) 

Furthermore, E g = —KN, where N is the number of electrons, and K = J2n=i l a n,d| 2 - 
Consequently, the starting Hamiltonian, H, becomes positive semidefinite 

H = P + U + E g , (5) 

except the additive constant E g . The transformation of the starting Hamiltonian based on 
the plaquette operators from Eq. (J2J) into Eq.(j3J) is possible only if Eq.® containing 19 
non-linear coupled equations allows solutions for the a UtP parameters. The different type of 
solutions of these non-linear equations are presented in the following. 

III. GROUND STATE SOLUTIONS 

We found that there are two types of solutions which satisfy the system of non-linear 
equations ©• These two types of solutions will be denoted by R — I and R = II. In /, 
for a* n d /a* n j — q n = q = real for all n, A- h „ reduces to one-site form A i (T = J2n=i a n,d^i+r„,o-) 
«4i,o- = d- 1)C + /i i(7 for all i. While in II, q n ^ q =real, and such a reduction of into A\ i(J 
is not possible. 

In the following we will analyze in details the ground state corresponding to the solutions 
/ and II: (A) We first determine the homogeneous phases at quarter filling (N = N\, 
being the number of lattice sites). (B) By decreasing N we find degenerated droplet, stripe 
and checkerboard ground states which we will present in detail, and (C) extensions of H 
are identified to lift the degeneracy leading to pure, non-degenerate stripe and checkerboard 
ground states. 

A. Ground states at quarter filling 

To find the ground state at quarter filling, a complementary unit cell operator— B i a is 
defined by 

{i iiCT ,%,} = 0, Vi,i>,<7'. (6) 
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For case II, B- ha = £0,0=0,1 J2 P =dj b n,pPl+r n ,^ [ as shown in Fig.fflb)], and defining n'(a, (3) = 
3 + a — (3 — 2aj3, and taking w = bi^/a^j, one finds 

b n j = -wa* n >,d, b n4 = +wa* n ,j . (7) 

While, for case /, B i<CT = T,ti=i b n ,dBi+r n ,a, and {A\. a , Bl a ,} = 0, i3 i)(7 = d- h(T - f i>a holds. 

Introducing D{ a {R) = (C{ a 8 R j + B\ a 5 RyII ), where C i)(T = {B\^ + vSl^) and v-, are 
arbitrary coefficients, the exact ground state at quarter filling becomes 

N A 

l* fl Ai/4> = IPU^)I°> , ^ = J,JJ (8) 

i=l 

where |0) is the bare vacuum. This is the ground state wave function, i.e., (P+U)\^ 9t R,i/4) = 
0. Here Ul^g^i/i) = since, for R = I the electrons with arbitrary spins are introduced 
on different sites, and for R = II all B- la in (JHJ) have the same fixed spin a, so the double 
occupancy for both p is excluded. Besides, based on (jHJ), further we have P\^f 9l R t i/4) = 0. 
Since the minimum eigenvalue of P + U is zero, thus 1^/3,^,1/4) is the ground state for non- 
zero, positive, although arbitrary U . The ground state for case I is degenerate and globally 
paramagnetic, with one electron on each site, i.e. the state is localized. While in case 77 
the ground state is non-degenerate, is a saturated ferromagnet, with 0, 1, or 2 electrons on 
any i site, e.g. the state is itinerant. 

B. Ground states below quarter filling 

1. Random droplet ground states 

Decreasing the number of D- ia operators in the product of Eq.fjHJ), the exact ground state 
below quarter filling can be written as well. In case J, for iV < N\, we obtain 

JV 

ivi/4>«>o) = n4i°>. ( 9 ) 

i=l 

where sites i can be arbitrarily chosen. Being empty sites present, the ground state is an 
itinerant paramagnetic phase. In case //, the ground state is reached only if touching B[ u 
operators defined on different sites [for example (5 19 o -, _B 2 5, CT ) in Fig. El a)] have the same 
spin in order to maintain the g ) = condition. These neighboring B[ a operators with 
fixed spin build up different blocks (droplets) Blj containing Nb^ particles. Two different 
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blocks have no common lattice sites, and their spin is non-correlated. If we denote by Nbi 
the number of blocks, then the ground state wave function becomes 

Nbi ( Nbi j \ 

lW>n>0>-ll P!Jl°). ( 10 ) 

where J2j ^Bi d — N, but otherwise, Nbij and the shape of the blocks remain arbitrary, as 
on the example shown in Fig. |2]b). This state is itinerant and represents a fully saturated 
ferromagnet up to iV > N\ — 8. In the limit N = N\ — 8 there is always a possibility for a 
new block to appear on the lattice, in a way that is not in contact with other blocks and as 
such it can have opposite spin. For example, in Fig. [2] a), the middle block -815,1 is isolated 
and has opposite spin compared to the big block surrounding it. Hence, for iV < N\ — 8 
the ground state is not anymore fully saturated. Further decreasing N, the ground state 
remains ferromagnetic^ until two disjoint blocks with the same number of sites but opposite 
spins can be constructed. This happens at iV? = N\ — 2L, where N\ = L x L. Thus, for 
N < Np, the ground state is globally paramagnetic. 



2. Degenerated stripe ground state solutions 

Decreasing N below N£, stripes emerge in the ground state in both Eqs. (|9I1U|) as vertical 
stripes (Fig. EJ) or diagonal stripes (Fig. 0J. The ground state wave function containing Nst 
vertical stripe lines is 

N St 

h«w*> = n( n ^Ui°>> ( n ) 

3=1 ieistj 

where Istj in case I (II) represents the stripe line j (plaquette stripe column j), and 
1st = J2j Ist,j- For example, from (fH|). the vertical stripes are obtained in case / (II) by a 
displacement along periodic vertical lines (vertically constructed plaquette columns) of C[ a 
(B\ ) operators. In Fig. El a) case I is shown with column stripes, while Fig. El b) depicts 
case II with plaquette columns. For diagonal stripes, the vertical displacement must be 
changed to diagonal one, as shown in Fig. HJ In the case II at dl = d2, see Fig. El the 
stripe structure turns into a checkerboard one. The stripe lines for case I are paramagnetic 
and insulating, while they are itinerant and ferromagnetic for case //. Different stripe lines 
which are not in contact have non-correlated spin. 
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C. Non-degenerated stripe and checkerboard ground states 

For N < Np, droplet [see Fig. |21 b)] and stripe solutions coexist as the ground state 
is degenerate. However, the droplet contributions can be eliminated in exact terms from 
the ground state by adding new Hamiltonian contributions to H. For example, let us 
consider Ha = — \W±\ Zae/ S t wnere ni = X^p^fcr If I'st contains all lattice sites from 
V\ lgt = Ylf =1 D ijtaj and 1st = I'st holds, then [E ie / St fk]v\, gt = T^JN + E ie / St fk] provides 
the minimum possible eigenvalue for Ha via Ha\^ g,R,i s ^l = — g^j^) . If however, 

I St 7^ 1st [for example, if j 2 moves to the in ji position, Fig. El a)], then \^ gt Rj St ) is not 
anymore an eigenstate of Ha- Consequently, \^ g ,Rj St ) becomes the unique, non-degenerate 
ground state of H + Ha- If we add to Ha the term H' A = \W 2 \ Hi^i st as well, the results 
remain unchanged. A Hamiltonian term of the form Ha is motivated in case of cuprates 
by low temperature tetragonal fluctuations^. The potential W in Ha = Ha + H' A can 
be generated by a periodic charge displacement or charge density wave (see Fig. 0), which 
is able to stabilize a stripe phase. Such a behavior has been already seen^S. in cuprates. 
Similarly, the checkerboard state also can be made stable^. A checkerboard is obtained, for 
example in Fig. E]b), if we take dl = d2 = 2. This is stabilized by a Hamiltonian term Ha, 
which has a set of lattice sites 1st, m such a way that it contains only next-nearest-neighbor 
sites on each second diagonal. 

Furthermore, the vertical stripes from Fig. El b), will be stabilized by a dimerization term 
of the form H B = E CT Eie/ E H^i^, E iiCT = J2 p =d,f{ e i >P Pi,a + e 2) pPi+ xl)Cr ), where I E contains 
only each second site of horizontal lattice rows^, see Fig. El For Hb described in Ref. |23|, the 
stable stripe phase occurs above the surface presented in Fig. [HI Modifying I E , other stripes 
can be addressed as well. The important effect of dimerization on stripes stabilization have 
been recently shown using numerical simulations^. 

IV. DISCUSSION OF THE RESULTS 

The obtained inhomogeneities are not connected in principle neither to certain special 
values of the parameters, nor to special properties of the homogeneous phases present at 
quarter filling. In the following we present several remarks to support this and in closing 
the obtained stripe formation process is summarized. 
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Firstly, for a given Hamiltonian, a given decomposition in positive semidefmite operators 
at a given filling factor is not uniqu e) 14 ' 25 ' 26 . The decomposition itself can be done in sev- 
eral different ways, leading to different A- X a operators, different matching conditions, hence 
different interdependencies between microscopic parameters obtained during the solution of 
the matching conditions (i.e. equations of the type 0), providing similar solutions, and 
similar characteristics in different regions of the parameter space. For example, as shown in 
Ref. instead of plaquettes, rhombi can also be used in defining A^ a , obtaining the same 
type of solutions, but in other regions of the parameter space. Other possibilities are to 
simply tilt the unit cell^, or to decompose in J2i A\ A, where Ai = S CT A,o-, as described in 
Ref.pS], etc. 

Secondly, the different solutions of the matching conditions all give homogeneous phase 
at quarter filling. Stripes will appear from each of these homogeneous phases with (hole) 
doping. This can be easily understood from the observation that the B i a (or B- ya ) operators 
which are characteristic of the stripe ground state, see eg, Eqs. (j6|7j) . can appear for arbi- 
trary A-^ a . Consequently, the obtained stripe formation process is weakly dependent on the 
properties of the homogeneous phase from which the stripes emerge. 

Thirdly, each homogenous phase (and hence each stripe structure originating from it) 
appears for different, well-defined conditions. For example, see also Ref. p^], stripes with 
ferromagnetic and itinerant stripe lines can be obtained from the itinerant ferromagnetic 
homogenous phase when the one particle part of the Hamiltonian in has such parame- 
ters that Hq will have a diagonalized partially filled lower flat band. One the other hand, 
insulating paramagnetic stripe lines are obtained in a different parameter space region, where 
localized homogenous phase can occur, etc. 

Fourthly, an additional repulsive Hubbard term acting in the <i-band will not change the 
obtained results. Consequently, the obtained ground states and inhomogeneities are valid 
not only in the 2D PAM, but also in a more general 2D two-band Hubbard model as well. 

Lastly, in our rigorous description we obtained a non-degenerate ground state exhibiting 
stripe structure in the following steps: we doped the homogeneous phase at quarter filling. 
This resulted in a degenerate ground state to appear, which contained both random droplets 
and stripe solutions. In the last step we lifted the degeneracy by adding a so-called stabi- 
lization term to the Hamiltonian. By doing this we eliminated the random droplets from 
the degenerated ground state, obtaining providing a non-degenerated stripe ground state 
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solution. 



V. CONCLUSIONS 

In conclusion, providing exact results for stripe, checkerboard vs. droplet interplay, we 
show how such intrinsic inhomogeneities appear in a 2D Hamiltonian as a non-degenerate 
ground states. For this, 1) a generalized PAM is used and cast in a positive semi definite 
form, 2) the ground states are explicitly constructed at and below quarter filling, and 3) the 
ground state degeneracy provided by random droplets is eliminated using extension terms 
representing, e.g., distortion lines, dimerization or density waves. The inhomogeneities were 
shown to exist in a broad concentration range below quarter filling and they are either 
(case I) paramagnetic and localized, or (case II) itinerant and ferromagnetic. In both of 
these cases the inter-stripe lines are insulating. As has been underlined, stripes can emerge 
from different homogeneous phases so are less sensitive to microscopic parameters of H. As 
marginal for the stripes, but important for PAM we derive an exact ferromagnetic ground 
state, as well, in 2D at quarter filling. An extra Hubbard interaction in the d band, Ud> 
leaves the above results unchanged. The presented positive semidefinite decomposition is 
not unique, can be effectuated in several ways, leading to similar results also in other regions 
of the parameter space. 
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FIG. 1: a)Unit cell (J), and, b) the B- 1S operator defined at an arbitrary site i, x T are the primitive 
vectors and n is the i independent notation of the sites in X\. For a) arrows indicate hopping and 
hybridization matrix elements (J=t, V). In b), the 6„ iP , p = d, f coefficients are shown together 
with a decomposition in horizontal components B- 1S = B^l+B^ X2 s of the complementary operator 
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FIG. 2: (a) Isolated plaquette with opposite spin at site i = 15, and (b) droplet solutions for case 
77. The blocks Bi yCr introducing the B\ a operators in the ground state wave- function of case II 
are presented as shaded plaquettes in b). 
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FIG. 3: Example of vertical stripe solutions for case I shown in (a) and II in (b). For a) (or b)) 
full circles (or shaded plaquettes) denote sites (or unit cells) whose 6\ a (or B\ a ) contribution is in 
the ground state wave-function Eq.(8). For example, site j2 in case I (unit cell (i2,i3,i5,i6) in case 
II). In both cases dl (d2) represents a measure of the stripe line width (inter stripe line distance) 
in x direction and |xi| units. 
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a) 



b) 



FIG. 4: Example of diagonal stripe solutions for case I shown in (a) and case II in (b). Notations 
are as in Fig.3. 
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FIG. 6: Bond alternation in x-direction to stabilize the stripes from Fig. 3 b), J x = t%, V a 




potential W created by Ha acting on the charge degrees of freedom in x direction stabilizing the 
stripes from Fig.3 a). 
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FIG. 8: (Color online) The Ef/\t^.\ surface in function of t = t^/t^ and the anisotropy parameter 
z = (l*y I ~~ I)/ I above which the R = II stripe ground-state from Fig.3.b emerge stabilized by 
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Hb m conditions from Ref.[23]. For hybridization V x = Vo = 0, and for other r, |T4|/|* 
holds. The inset shows the required tf~ alternation introduced by Hb- For Hb = 0, FM is present 
on the plotted surface at 1/4 filling. 
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